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Abstract 

Single and pair top-quark production in proton-lead (p-Pb) and lead-lead (Pb-Pb) collisions at the CERN Large 
Hadron Collider (LHC) and future circular collider (FCC) energies, are studied with next-to-leading-order pertur¬ 
bative QCD calculations including nuclear parton distribution functions. At the LHC, the pair-production cross 
sections amount to utt = 3.4 ph in Pb-Pb at = 5.5 TeV, and = 60 nb in p-Pb at y/s^ = 8.8 TeV. At the 
FCC energies of = 39 and 63 TeV, the same cross sections are factors of 90 and 55 times larger respectively. 
In the leptonic final-state ti —>■ W^b W~b b b££ vv, after typical acceptance and efficiency cuts, one expects 
about 90 and 300 top-quarks per nominal LHC-year and 4.7-10"^ and 10® per FCC-year in Pb-Pb and p-Pb col¬ 
lisions respectively. The total ti cross sections, dominated by gluon fusion processes, are enhanced hy 3-8% in 
nuclear compared to p-p collisions due to an overall net gluon antishadowing, although different regions of their 
differential distrihutions are depleted due to shadowing or EMC-effect corrections. The rapidity distributions of 
the decay leptons in ti processes can be used to reduce the uncertainty on the Ph gluon density at high virtual¬ 
ities by up to 30% at the LHC (full heavy-ion programme), and hy 70% per FCC-year. The cross sections for 
single-top production in electroweak processes are also computed, yielding about a factor of 30 smaller number of 
measurable top-quarks after cuts, per system and per year. 


1 Introduction 


The multi-TeV energies available at the CERN Large Hadron Collider (LHC) have opened up the possibility to 
measure, for the first time in heavy-ion collisions, various large-mass elementary particles. After the first observa¬ 
tions of the W 00 and Z @0 bosons, as well as bottom-quark (6-jets) 0 ^ there remains only three Standard 
Model (SM) elementary particles to be directly measured in nucleus-nucleus collisions: the r lepton, the Higgs 
boson, and the top quark. Whereas the r measurement should be straightforward, that of the Higgs boson is be¬ 
yond the LHC reach as it requires much larger cross sections and/or luminosities |j^, such as those reachable at the 
proposed future circular collider (LCC) Q with about seven times larger center-of-mass energies than at the LHC. 
The study presented here shows, for the first time, that the top-quark -the heaviest elementary particle known- will 
be produced (singly or in pairs) in sufficiently large numbers to be observed in lead-lead (Pb-Pb) and proton-lead 
(p-Pb) collisions at the LHC and LCC. 

Since the width of the top-quark (Ft « 2 GeV) is much larger than the parton-to-hadron transition scale given 
by « 0.2 GeV, the top-quark is the only coloured particle that decays before its hadronization. Its short 

lifetime, tq = /i/Ft « 0.1 fm/c, implies that the top decays -into a t W b final-state with a nearly 100% branch¬ 
ing ratio 0- mostl}j^ within the strongly-interacting medium, such as the quark-gluon plasma (QGP), formed in 
nuclear collisions. The large top-quark mass (rrit « 173 GeV) provides a hard scale for high-accuracy perturbative 

' The typical transverse momentum of the produced top quark is usually smaller than its mass, px < rrit, and the Lorentz-hoost factor is 
7 Ri cosh [yt), where yt is the t-quark rapidity. At the LHC (\yt\ < 3) the Lorentz-dilated mean decay time is r = 7 T 0 ~ 0.1 — 1 fm / c, and at 
the FCC (Is/t I < 5), r = 7 T 0 ~ 0.1 — 7.5 fm/c; to be compared with the typical QGP formation time of 1 fm/c, and lifetime of 10 — 20 fm/c. 
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calculations of its quantum-chromodynamics (QCD) and electroweak production cross sections (next-to-next-to- 
leading-order, or NNLO, is the current theoretical state-of-the-art sini)- At hadron colliders, top quarks are 
produced either in pairs, dominantly through the strong interaction, or singly through the weak interaction. At the 
energies considered here, the dominant production channels, as obtained at NLO accuracy with the MCFM code p^ , 
are: (i) gluon-gluon fusion, g g ^ tt + X, contributing by 80-95% to the total pair production (the remaining 5- 
20% issuing from quark-antiquark annihilation), (ii) f-channel single-top electroweak production qb ^ q't + X 
(the s-channel process, decreasing with energy, amounts to 5-1.5% of the total single-f cross section), and (iii) 
associated top plus kk-boson, gb ^ Wt + X, production (increasing with energy, it amounts to 25-50% of the 
f-channel process). 

The theoretical motivations for a dedicated experimental measurement of the top-quark in heavy-ion collisions 
are varied and include, at least, the following studies: 


(i) Constraints on nuclear parton distribution functions (nPDFs). The tt cross sections in proton-proton (p-p) colli¬ 
sions can be used to constrain the proton PDFs GD- In the heavy-ion case, top-pair production probes the nu¬ 
clear gluon density in an unexplored kinematic regime around Bjorken-a: values, x « 2 rut /y/s^ ~ 5T0“^- 
0.05, and virtualities « 3 • 10"*^ GeV^, a region characterized by net positive, albeit small, anti¬ 

shadowing corrections. In addition, at the FCC, the &-quark nPDF (in single-top production), and even the 
top-quark nPDF itself, are generated dynamically by the constituent gluons and become necessary ingredients 
of the theoretical cross section calculations. 


(ii) Heavy-quark energy loss dynamics. The top quark can radiate gluons before its W b decay which, given its 
very-short lifetime, occurs mostly inside the QGP Medium-induced gluon radiation off light-quarks and glu¬ 
ons, leading to “jet quenching” p^ , results in a factor of two reduction of jet yields in Pb-Pb compared to 
' = 2.76 TeV|[l5l|16|. Although solid theoretical expectations for heavy-quark radiation 


p-p collisions at 


/ c — 


predict a reduced amount of gluonstrSilung at small angles due to the “dead cone” effect HI). the experi¬ 
mental data somehow unexpectedly shows the same amount of suppression for jets from light-flavours and 
6-quarks 0. The relative role of elastic and radiative scatterings on the energy loss of heavy-quarks is an 
open issue in the field 118 The detailed study of top-quark production in heavy-ion collisions would 
therefore provide novel interesting insights on the mechanisms of parton energy loss. In addition, the study 
of boosted top-pairs (with transverse momenta above px ~ 1 TeV) traversing the QGP as a colour-singlet 
object for a fraction of their time, will allow one to probe the medium opacity at different space-time scales. 


(iii) Colour reconnection in the QGP. The top mass, featuring the strongest coupling to the Higgs field, is a funda¬ 
mental SM parameter with far-reaching implications including the stability of the electroweak vacuum p0| . 
Currently, the dominant rrit systematic uncertainty is of theoretical nature and connected to the modeling of 
the colour connection and QCD interferences between the tt production and decay stages, and among the 
hadronic decay products. Indeed, the colour-flow (through gluon exchanges and/or non-perturbative string 
overlaps) between the t and t quarks, their decayed 6-quarks, and the underlying event from multi-parton 
interactions and beam-remnants sutTounding the initial hard scattering pT) , results in uncertainties on the re¬ 
constructed rrit of a few hundred MeV. The amount of top quark interactions with the colour fields stretched 
among many partons involved in nuclear collisions will be obviously enhanced compared to more elementary 
systems. Thus, the reconstruction of the top-quark mass in the QGP (assuming its feasibility is not jeopar¬ 
dized by the large 6-quark energy loss already observed in the data), or in proton-nucleus interactions, would 
provide interesting insights in non-perturbative QCD effects on a crucial SM parameter. 


In this paper we mostly focus on nPDF constraints through top-pair production in p-Pb and Pb-Pb collisions. 
We also provide the expected px reach of top quark spectra at various energies to indicate where boosted final-states 
can be measured for energy loss studies. The paper is organized as follows. In Sectionj^the theoretical setup used is 
outlined, which is then used to compute the NLO cross sections at the LHC and FCC, and associated yields expected 
after typical acceptance and efficiency cuts, for top-pair and single-top production, presented in Section]^ Sectionj^ 
quantifies the impact on the nuclear PDFs provided by the measurement of the rapidity distributions of the decay 
leptons from top-quark pairs produced at the LHC and FCC, using a Hessian PDF reweighting technique p^[23| . 
The main conclusions of the work are summarized in Section|5] 
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2 Theoretical setup 


The top-pair and single-top cross sections are computed at NLO accuracy with the MCFM code (version 6.7) 
using the NLO CTIO proton PDFs (including its 52 eigenvector sets) p4| corrected for nuclear effects (shadowing, 
antishadowing and EMC) | [25| through the EPS09 nPDEs (including its 30 error sets) for the Pb ion As our 
main purpose is to provide estimates for the feasibility of different top-quark measurements in nuclear collisions, we 
do not discuss here the (subleading) sensitivity of the cross sections to different sets of proton PDEs (nor associated 
variations of the strong coupling as). Also, while there are other nuclear PDE sets available p7|j29) , we only 
employ EPS09 here as it is the only nPDE set that is consistent with the dijet measurements in p-Pb collisions at 
the LHC 130 ^ (the data would probably agree also with the latest nPDEs by nCTEQ p2) , but these sets are not 
available at the time of writing this article). We run the following MCFM processes: i4i for tt production, lei, i66 for 
single-(anti)top in the f-channel, ivi, i7 6 in the s-channel, and lei, i86 for associated t W production. We note that, 
at NLO, the theoretical processes defining t W production partially overlap with those contributing to top-quark 


pair production 1331. This is accounted for in our MCFM t W cross sections calculations by vetoing the additional 
emission of a 5-jet. The code also properly accounts for the different isospin (u- and d-quark) content of the Pb 
nucleus, which has a small impact on the electroweak single-top processes. 

All numerical results have been obtained using the latest SM parameters for particle masses, widths and cou¬ 
plings 1^, and hxing the default renormalization and factorization scales at p = /ij? = = rrit for tt and t-,s- 

channel single-top, and at p = ~ t W processes. The NLO calculations used 

here reproduce well the cross sections experimentally measured at the LHC in p-p collisions at = 7, 8 TeV for 
tt 134 ■ ^, f-channel single-top 1 46p9 -5^, and associated tW \ 53 - ^ production. Incorporation of next-to-NLO 
corrections 0 would increase the theoretical cross sections, by about 10%, i.e. the so-called AT-factor amounts to 


K = cr,, 


)A 


1.10, and further improve the data-theory agreement. The computed nucleon-nucleon cross 


sections are scaled by the Pb mass number {A = 208) to obtain the p-Pb cross sections, and by = 43 264 in 
the Pb-Pb case, as expected for hard scattering processes in nuclear collisions. The uncertainties of the theoretical 
cross sections are obtained from the values computed using the eigenvector sets of, hrst, the CTIO and, then, EPS09 
PDEs and adding them in quadrature, as well as by independently varying the renormalization and factorization 
pLp and pp scales within a factor of two (for the central CTIO and EPS09 set). The PDE and scale uncertainties 
amount each to 3-10% for tt and 3-6% for single top. The scale uncertainties have no impact on the nPDE con¬ 
straints derived below given that the scale dependence of the nuclear modihcations in PDEs around /r = mt is very 
mild (see e.g. Eig. 1 in Ref. | [5^ ) and that differences between proton and nuclear PDEs are obtained via ratios of 
(p-Pb,Pb-Pb)/(p-p) cross sections at the same colliding energy, where those mostly cancel out. 


3 Top-pair and single-top cross sections and yields 


Table collects the total cross sections, and associated scale and PDE uncertainties, for top-pair and single-top 
production in p-Pb and Pb-Pb collisions at LHC and ECC energies, obtained as described in the previous Section. In 
the case of ti and t W production, a net EPS09 gluon antishadowing in collisions with Pb ions results in an increase 
of the total production cross sections by about 2-8% compared to those (A-scaled) obtained for p-p collisions using 
the proton CTIO PDE. Eor t- and s-channel single-top cross sections, the overall nuclear modifications are quite 
insignificant, ±2% depending on the energy. Eigure shows the total top-pair and single-top cross sections as a 
function of collision energy for p-p, p-Pb and Pb-Pb collisions in the range of c.m. energies in the nucleon-nucleon 
system of y^sj^ « 1-100 TeV. The single-top curves are obtained adding the f-,s-channel and tW cross sections 
listed in Table [T] In general, top-quark pair production is a factor of two (four) larger than the sum of single top 
processes at the LHC (ECC). The nominal LHC and ECC energies for p-Pb and Pb-Pb collisions are indicated by 
dashed boxes in the plot. Going from LHC to ECC, the total cross sections increase by signihcant factors, x (55-90) 
for tt and x (30-40) for single-top. 

The impact of nuclear PDE modihcations on the yields for a given hard process is usually quantihed through the 
nuclear modihcation factor Rj^j^ given by the ratio of cross sections in nuclear over proton-proton collisions scaled 
by A or A^. The theoretical (l/t.t) factors as a function of the rapidity of the produced top and antitop quarks 
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Table 1: Inclusive cross sections for top-pair and single-top (f-channel, s-channel, and tW) production in Pb-Pb 
and p-Pb collisions at LHC and FCC energies, obtained at NLO accuracy with MCFM. The first uncertainty is due 
to theoretical scale variations, and the second one to the CTIO and EPS09 PDF errors added in quadrature. 


System / Process: 

(MCFM process) 

top pair (tt) 

(l4l) 

single-top (f-channel) 

(l61,166) 

single-top (s-channel) 

(l71,176) 

single-top {t W) 

(l81,186) 

Pb-Pb = 5.5 TeV 

p-Pb ys;^=8.8TeV 

3.40±0.42±0.37 /ib 

58.8±7.1±3.8nb 

L61±0.05±0.08 pb 

2L1±0.63±0.63 nb 

110±4±6nb 

L09±0.03±0.04nb 

313±13±41 nb 

5.26±0.21±0.37nb 

Pb-Pb yv; = 39 TeV 
p-Pb = 63 TeV 

302±33±12 pb 

3.20±0.35±0.10/ib 

54.6±L6±2.2/ib 

518±16±17nb 

L31±0.05±0.08 pb 

10.9±0.5±0.5 nb 

24.2±L6±L3 pb 

246±24±11 nb 



Figure 1: Total cross sections for top pair and single-top (sum of f-,s-channels plus t W processes) production in 
Pb-Pb p-Pb and p-p collisions as a function of c.m. energy. The dashed boxes indicate the nominal nucleon-nucleon 
c.m. energies, = 5.5, 8.8, 39, 63 TeV, of the heavy-ion runs at the LHC and FCC. 


4 























are shown in Fig. [^for tt at LHC (top panels) and FCC (bottom panels) energies. The central curves indicate the 
result obtained with the central EPS09 set and the grey bands show the corresponding nPDF uncertainties. All the 
results presented in those, and the following, plots are given in the center-of-mass frame of the colliding species. 
In general, the distributions reveal very similar trends for ti and single-top (not shown here) processes, 

although the f-channel and s-channel processes have a factor of two smaller uncertainties, as expected, given that 
single-top production is dominated by quark-induced processes whose densities in the nucleus are better known than 
the gluon ones which produce most of the top pairs. At the LHC, for both (single and pair) production mechanisms, 
the nPDF effects increase the average top-quark distributions at central rapidities by about 10% (antishadowing) 
while they deplete them by 20% at backward rapidities (also in the forward direction in Pb-Pb) due to the so-called 
EMC effect at large-a: p5| (see also Eig. [plater). At the ECC, the higher collision energies as well as the larger 
kinematical coverage assumed for the detectors at this future facility give access to smaller momentum fractions x 
where EPS09 predicts moderate shadowing even at high virtualities Q ^ mt- This leads to additional suppression 
at forward rapidities (and also in the backward direction in Pb-Pb). 



Eigure 2: Theoretical nuclear modification factor as a function of rapidity of the individual top and antitop quarks in 
tt production in Pb-Pb (left panels) and p-Pb (right panels) at = 5.5, 8.8 TeV (LHC, top panels) and 39, 63 TeV 
(ECC, bottom panels), computed at NLO accuracy with MCFM. The central curve indicates the result obtained with 
central EPS09 parametrization and the grey band the corresponding nPDE uncertainty. 

The cross sections listed in Table[T]and plotted in Eig.[T]are total inclusive ones and do not include the t ^ W b 
decays, nor any experimental acceptance/analysis requirements on the final-state particles. The determination of 
the expected yields at the LHC and ECC requires accounting for top and PF-boson decays plus acceptance and 
reconstruction efficiency losses. The W leptonic branching fractions, W —>■ £^i/( (with i = e, y, r), amount to 
1/9 for each lepton flavour, the other 2/3 being due to W dijet (quark-antiquark) decays. In this work we will only 
consider leptonic W decays characterized by a final-state with an isolated electron or muon plus missing transverse 
energy (1^^) from the neutrino, because the W dijet-decays are much more difficult to reconstruct in the large 
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background of heavy-ion collisions, and also potentially subject to parton energy loss effects (although they could 
be certainly tried in the “cleaner” p-Pb environment). For the top-pairs measurement, ti —> W^h W~b b b£i vv, 
the combination of electron and muon decays for both FF-bosons (ee, /i/r, e/r, /re) reduces the total cross sections 
by a facto][^of = 4/9^ « 1/20. For the single-top case, we will only consider the yields for associated tW 
production which shares a characteristic final-state signature, tW W bW b££ vv, very similar to that of 
top-pair production. Indeed, the experimental observation of f-channel (let alone the much more suppressed s- 
channel) single-top, with one less charged lepton and neutrino, is much more challenging on top of the expected 
large W, Z+jets background (in Pb-Pb, at least, although it should be feasible in p-Pb collisions). 


In order to compute the expected number of top-quarks measurable at the LHC and FCC, we include in the 
MCFM generator-level calculations the typical analysis and fiducial requirements for 6-jets, charged leptons, and 
missing transverse energy from the unidentified neutrinos, used in similar p-p measurements at the LHC | L5|38|54 |. 
Although some of these p-p experimental requirements may seem optimistic for the more complex environment 
encountered in heavy-ion collisions, they are validated by future experimental projections of the CMS collabora¬ 
tion In the case of FCC, we extend the pseudorapidity coverage for charged-lepton tracking and 6-jet secondary 
vertexing from the LHC range of jpl = 2.5, to [t/I = 5. The details of all selection criteria are given in Table|^ We 
reconstruct the 6-jets with the anti-A:.j, jet clustering algorithm | [5^ with distance parameter R = 0.5, and we require 
the high-pT charged lepton to be separated from the closest 6-jet within an (?/, (j)) isolation radius of = 0.3. 


Table 2: List of analysis cuts on single px and p for 6-jets and isolated leptons £ = e^, and on the neutrinos 


, typically employed in top-pair (tt —>• 6 h££vv) |35 38 and single-top plus IL-boson (/FF —>• b££vv) |54 


measurements in fully-leptonic final-states in p-p collisions at the LHC, applied in our generator-level studies. 


Analysis cuts 

6-jets (anti-k^ algorithm with R = 0.5): px > 30 GeV; \r]\ < 2.5 (LHC), 5 (FCC) 
charged leptons £ (R,,.„i = 0.3): px > 20 GeV; \r]\ < 2.5 (LHC), 5 (FCC) 
neutrinos: > 40 GeV 


The combination of all analysis cuts listed in Table [^results in total acceptances of order « 40% 

(50%) for tt and fFF measurements at the LHC (FCC). In addition, one has to account for experimental 6- 
jet tagging efficiencies, which we conservatively take of the order of 50% as determined in Pb-Pb collisions at 
^/s = 2.76 TeV ||^. For single-top, this results in an extra « 0.5 reduction of the measured yields. In the ti case, 
in order to tag the event as such, one usually only requires a single 6-jet (out of the two produced) to be identified, 
and thus the associated efficiency is larger: « 1 — (1 — 0.5)^ = 0.75. The combination of acceptance, analysis 

requirements, and efficiency losses results in an overall efficiency factor of x « 30% (40%) for the final tt 
yield^at the LHC (FCC). The same factor for t W production is A^.^^ x « 20% (25%) at the LHC (FCC). Pos¬ 
sible backgrounds, mostly from FF, Z+jets, WZ, and ZZ production sharing similar final-state signatures as both 
top-quark production processes, can be minimized by applying dedicated jet-veto requirements and/or extra criteria 
on the invariant masses of the two high-px leptons, e.g. away from the Z boson peak {\m^ — mu\ > 15 GeV). 
We do not directly compute the impact of such backgrounds here as the applied analysis cuts already reduce them 
to a manageable level according to the existing p-p measurements. In particular, the application of mu cuts would 
reduce the visible yields by an extra 10% which is, however, compensated by the fact that our NLO calculations 
need to be scaled by about the same amount to match the current p-p data (and NNLO predictions). We note also 
that, despite the larger hadronic backgrounds in nuclear collisions, the instantaneous luminosities in the heavy-ion 
operation mode at the LHC (and FCC) result in a very small event pileup, at variance with the p-p case, and make 
the top-quark measurements accessible without the complications from overlapping nuclear collisions occurring 
simultaneously in the same bunch crossing. 

^Including also e=*= and from leptonic tau-decays in the tt -+ er, tt + final-states would decrease the con'esponding branching 
ratio only to Pi 1/16. 

^We note that, in the Pb-Pb case, parton energy loss effects which can bring the 6-jet below the px threshold criterion (px > 30 GeV) and 
result in an additional inefficiency to tag the tt event, are unlikely to affect both 6-jets at the same time. Indeed, for simple geometrical reasons 
if one top-quark is produced and decays through the denser region, the other one emitted back-to-back will go through a thinner medium layer 
and its associated 6-jet will be tagged with our considered probability. 
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The expected total number of top-quarks (adding the t and i values) produced in Pb-Pb and p-Pb in one year at 
the nominal luminosities for each colliding system, obtained via N = a ■ ■ A - e, are listed in Table The 

number of visible single tops in tW processes is smaller by a factor of ^30 compared to those from ti production, 
due to a combination of causes: lower cross sections, smaller reconstruction efficiencies, and only one top-quark per 
event. At the LHC, we expect about 100 and 300 top-quarks measurable in the fully leptonic decays from tt-pairs 
and t W processes in Pb-Pb and p-Pb respectively. For comparison, the CDF and DO experiments reconstructed less 
than 100 top-quarks (in all decays channels) during the full Run-1 operation at Tevatron. At the end of the LHC 
heavy-ion programme, with « 10 nb~^ (1 pb“^) integrated in Pb-Pb (p-Pb), about 2.5 thousand (fully-leptonic) 
(anti)top-quarks should have been measured individually by the CMS and ATLAS experiments. The corresponding 
visible yields at the FCC are about 300 times larger, reaching 5 x 10"* and 10® top-quarks per year in Pb-Pb and p-Pb 
collisions respectively. 

Table 3: Expected number of top-rantitop quarks per run, after typical acceptance and efficiency losses, for top- 
pair and t W single-top measurements in fully-leptonic final-states in p-Pb and Pb-Pb collisions at LHC and FCC 
energies for the nominal per-year luminosities quoted. 


System 


r 

^int 

Number of top-rantitop quarks 

ti ^ bhtivv 

Number of top-rantitop quarks 

tW ^ hUvv 

Pb-Pb 

5.5 TeV 

1 nb-® 

90 

3 

p-Pb 

8.8 TeV 

0.2pb-i 

300 

10 

Pb-Pb 

39. TeV 

5nb-® 

47 000 

1300 

p-Pb 

63. TeV 

1 pb-® 

100000 

2 600 


In order to provide an idea of the top-quark px reach accessible in the different measurements listed in Table 
we have also computed the expected t, t transverse-momentum distributions in Pb-Pb and p-Pb collisions, after all 
acceptance and efficiency criteria applied (Fig.[^. The maximum top-quark px experimentally measurable per LHC 
(FCC) year will be around px 300 (1500) GeV for Pb-Pb, andpx 500 (1800) GeV for p-Pb. Given the limited 
LHC top-quark statistics, the study of boosted-tops will be thus only accessible at the future circular collider. 


4 Constraints on nuclear PDFs from tt production 


As aforementioned, 80-95% of the total pair production at LHC-FCC comes from gluon-gluon fusion processes 
and, thus, ti cross sections can be used to constrain the relatively badly-known gluon densities in the Pb nucleus. 
In this section we quantify the impact that top-quark measurements at the LHC and FCC would have on bet¬ 
ter constraining the nuclear PDFs through the so-called Hessian PDF reweighting technique p2l[2^ . Such PDF 
reweighting procedure is based on the fact the error sets {/}J defined in Hessian PDF-fits correspond to a certain 
increment ( Ax^ = 50 forEPS09) of the global “goodness-of-fit” function whose minimum Xo is achieved 
with the central set {/}o- The error sets thereby constitute a parametrization of the original x^ function which can 
be taken advantage of in order to determine the associated PDE uncertainty after adding new experimental datasets. 
More precisely, the Hessian method for determining the PDE errors writes the response of the original 
function to fit-parameter variations 6ai as 

X^{a} - Xo ~ X] zl, ( 1 ) 

i,j k 

where the Hessian matrix Hij = {l/2)d^x^/(daidaj) is diagonalized in the last step. The PDE error sets {/}^ are 
then defined by Zi{{f}^) = ±^/A^Sik■ The impact of including new experimental data can now be computed 
by considering a function 

xLw ^ E +E ’ (2) 

k i,j 
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Figure 3: Expected top-quark px distributions, /dptf, in Pb-Pb (left panels) and p-Pb (right panels) in the 

fully-leptonic decay modes at = 5.5, 8.8 TeV (LHC, top panels) and 39,63 TeV (FCC, bottom panels) after 
acceptance and efficiency cuts. The curves are a fit to the underlying MCFM distribution. The markers indicate 
pseudodata corresponding to the luminosities listed in Table 
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where Di denote the ith new data-point and C is the covariance matrix that encodes the experimental uncertainties. 
The corresponding theoretical values are denoted by Ti and they depend on the PDFs. To first approximation, the z 
dependence of each Ti is given by 

T, (z) ^ r, (0) + ^ ~ ^ Wk. (3) 

k 

where Wk = ZkjTi (0) is the theory value co mp uted with the central set {/}o, and Ti{z^) is the theory 
value evaluated with the error set {f}t- Inserting Eq. ([^ into Eq. ([^ and requiring dXnev^/dwk = 0 for all fc, one 
finds the condition = -J2i B^iO^i for the minimum of xLw’ where Bkn = J2i,j DikC~^^Djri + Ax^Skn , 
DikC~j^ [Tj (0) - Dj], and Dik = - Ti {z^)) /2. The new theory values are obtained 

directly from Eq. ([^ with Wk = and the corresponding new central set of PDEs {/"®'^}o by replacing Ti 

with the PDEs, Ti —> Q^). To find the new PDE error sets we rewrite Eq. (j^ as 

xLw - xLw I = X] Bij Swj=J2'^l^ (4) 

i,j k 

where 5wi = Wi — and in the last step the matrix B is being diagonalized. The new PDE error sets 

can be then defined exactly as earlier by The procedure sketched above was proven 

in Ref. [ [2^ to be consistent with the Bayesian reweighting method introduced originally in Ref. | [60) and further 
confirmed more recently in Ref. 0- 


To mimic a realistic experimental situation, we generate sets of pseudodata for nuclear-modification factors 
and -i.e. for the ratios of cross sections obtained with EPS09 nPDEs over those obtained with the 

CTIO proton PDEs- corresponding to the LHC and ECC scenarios discussed in the previous Section. As noted 
earlier, the total tt production cross section is expected to undergo only a mild increase due to nuclear effects in 
the nPDEs, and it may be challenging to resolve such an effect from overall normalization uncertainties. Thus, the 
total tt cross sections are not expected to have as large impact as they have on the absolute free proton PDEs 1131. 
In order to get better constraints on the nPDEs differential cross sections are thus needed. Eor this purpose, we 
use the distributions of leptonic top-decay products, which are unaffected by final-state interactions in the strongly- 
interacting matter produced in nuclear collisions. Here, we concentrate on the measurement of tt pairs via the 
tt —t bb + + vv decay channel, binned in the charged lepton rapidity dNi/dyi with £ = e, y. The baseline 

for the pseudodata is taken from the nuclear modification factors computed with the central set of EPS09 in the 
previous Section. The expected number of events N (Ay^) in each rapidity bin Ayi are computed by 


N (Ayi) = iVtotai X 


a{yi G Ay) 

^total 


(5) 


where Atotai is the total number of events expected after acceptance and efficiency losses for each systenj^ listed 
in Table a{yi G Ay) is the cross section within rapidity bin Ayi, and Ctotai is the total cross section within 
the acceptance. The statistical uncertainty is then taken to be i/lJWjAyij to which we add in 

quadrature a constant ±5% systematic error (<5^^’’* = 0.05 x such that the total uncorrelated error is 

^uncorr _ ^ ^jstat ^2 _|_ 'pjjg overall normalization error is taken to be 5% (5™''™ = 0.05 x j^EPsog^ 

Systematic uncertainties of this order are realistic as the corresponding p-p measurements at the LHC 138 ^ 
have already reached a better precision. In addition, partial cancellations of systematic uncertainties common to 
p-p, p-Pb, and Pb-Pb measurements are expected in a careful experimental determination of and i?pi,pb- 
The statistical precision of the p-p baseline data, taken at slightly different center-of-mass energies, is also much 
better than that of the p-Pb and Pb-Pb measurements, and the additional theory uncertainty for energy-dependent 
corrections of the reference cross sections (e.g. to go from ^/s = 8 TeV in p-p to = 8.8 TeV in p-Pb at 
the LHC) are small 1^ . Each pseudodata point Di is then computed from the baseline values T^psog fj-om 
the estimated uncorrelated and normalization errors by Di = _|_ juncorr^,^ _|_ jnorm^norm^ ^ where r; and 

^norm random numbers from a Gaussian distribution of variance one centered around zero. The elements of the 
covariance matrix C are computed as Cy = [5“’ 


i^uncorr runcorr r 

Uj Uij 


jTnorm fnorm 


|63|. 


^For the LHC, we consider the total luminosity to be accumulated during the full heavy-ion programme, which is a factor of 10 (5) higher 
than the nominal per-year luminosities quoted for Pb-Pb (p-Pb). 
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Figure 4; Impact on the nuclear glue of tt (pseudo)data for the full LHC heavy-ion programme. Left; Nuclear 
modification factor for tt decay-leptons as a function of rapidity, Rp^pi, (ye), obtained for Pb-Pb at = 5.5 TeV 
(top) and p-Pb at = 8.8 TeV (bottom) compared to predictions computed with EPS09; current nPDF set 
(region enclosed by the red dotted lines) and after pseudodata-reweighting (blue curve plus grey band). Right; 
Ratio of nuclear-over-proton gluon densities, R^ evaluated at Q = mt, for the original EPS09 (band enclosed by 
red dotted lines) and for the reweighted EPS09 (blue curve with grey band) for Pb-Pb (top) and p-Pb (bottom). 


In the original EPS09 analysis p^ , the inclusive pion data measured at RHIC | [64) was given an additional 
weight factor of 20 in the x^-function in order to enhance the constraints on the badly known gluon densities from 
nuclear deep-inelastic and fixed-target Drell-Yan data alone. As both the pion data and the top quark production 
considered here, are mostly sensitive to nuclear gluon PDEs, we rescale the covariance matrix equally by C —C/20 
when performing the reweighting. This compensates for the large weight given for the RHIC pion data and should 
lead to a more realistic estimate of the impact that the top-quark measurements would have if directly included into 
the EPS09 fit. After finding the new theory values T”®"' through the reweighting procedure, the optimum overall 
shift (originating from the allowed uncertainty in normalization) is found by solving the multiplicative factor / from 
the contribution of the new data (see e.g. Ref J^); 


E 


[T— - D,] C-/ [t;— - D,] = min ^ 




-D,- fSf' 


1 2 


+ f 


( 6 ) 


In the results presented below, the resulting shift has been applied on the data points. 

The results of the nPDE reweighting procedure are presented first for LHC energies in Eigurej^for Pb-Pb (top) 
and p-Pb (bottom) collisions. The nuclear modifications for the decay leptons are somewhat less pronounced in 
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comparison to the top distributions themselves (Fig. ^ due to the smoothing brought about by the additional phase- 
space integrations related to the top-quark decays. The estimated statistical errors are generally of the order of 
10% and are foreseen to dominate over the predicted systematic uncertainties. These pseudodata probe the nuclear 
PDFs predominantly in the range 0.03 ^ x < 0.3 at ^/s = 5.5 TeV, and 0.02 < a; < 0.2 at y/s = 8.8 TeV, as 
inferred from our MCFM calculations. In both cases the pseudodata are found to have only a moderate impact on 
the EPS09 gluon densitj0 This is predominantly due to the rather low foreseen statistics and the rapidity interval 
covered, which makes especially somewhat flat within the acceptance. Consequently, even the overall 

normalization alone can mimic the effects of nuclear PDFs thereby reducing the obtainable constraints. The new 
error bands in Figs.|^are both around 10% narrower than the original EPS09 ones. Combining the p-Pb and Pb-Pb 
measurements and assuming independent data samples available from both the CMS and ATLAS collaborations, the 
total impact of tt production on the large- and mid-x nuclear gluons could reach 30% (with the full LHC luminosity). 
Such a modest improvement will most likely be overpowered by the constraints offered by the inclusive jet and dijet 
data from the LHC p-Pb run(s) pT] . 




10 


10 ' 


10 


10 




10 


10 '‘ 


10 


10 




Eigure 5: As Eigure|^but for the case of Pb-Pb (top) and p-Pb (bottom) per ECC-year pseudodata. 


The results of repeating the reweighting procedure with the ECC pseudodata are shown in Eig.|^for Pb-Pb (top) 
and p-Pb (bottom) collisions respectively. Although the foreseen ECC per-year integrated luminosities are similar to 
those expected for the full LHC heavy-ion programme, the production cross sections being a factor x (55-90) above 
the LHC expectations (Table signihcantly increase the expected top-quark yields thereby reducing the statistical 
uncertainties. Indeed, in our ECC scenario the systematic uncertainties dominate. The reduced uncertainties as well 

^As the dependence of Itg^(x, Q^) is rather mild for > 10 GeV^ |56| , the plots at are representative for most practical 

applications. 
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as the wider kinematic reach at the FCC make the impact of these pseudodata on EPS09 clearly larger than that 
expected at the LHC. The constraints also reach lower values of x, the dominant x region being 5 x 10“"* x < 
3 X 10“^ at y/s = 39 TeV, and 2 x 10“^ ^ 2; < 2 x 10“^ at yfs = 63 TeV. The addition of the nuclear ti results 
shown in Fig.|^would allow one to decrease the gluon density uncertainty by up to 50% in some regions. Unlike in 
p-Pb collisions, in the Pb-Pb case the nuclear effects coming from small-x (shadowing) and large-x (EMC effect) 
cannot be distinguished since they are essentially multiplied at large \yi\. In this case, the new constraints tend 
to affect more the part that is originally less constrained, the high-x gluons. That is, while the probed region in 
X is very similar in our p-Pb and Pb-Pb scenarios, the p-Pb collisions are foreseen to provide more varied nPDF 
constraints. Combination of p-Pb and Pb-Pb data, plus assumption of two independent experiments measuring the 
spectra, would result in an overall reduction of up to 70% with just one year of integrated luminosity. 


5 Conclusions 


The study presented here has shown, for the first time, that top quarks produced in pairs via (mostly) gluon-gluon 
fusion, or singly in electroweak processes, are clearly observable in p-Pb and Pb-Pb collisions at the energies of 
the CERN LHC and future circular collider (FCC). The corresponding cross sections have been computed at NLO 
accuracy with the MCFM code including the CTIO free proton PDFs and nuclear modifications parametrized with 
the EPS09 nPDE At the LHC, the pair-production cross sections are = 3.4 pb for Pb-Pb at = 5.5, and 
a-tt = 50 nb for p-Pb at = 8.8 TeV. At the FCC energies of = 39, 63 TeV, the same cross sections are 
factors of 90 and 55 times larger respectively. The total tt cross sections are enhanced by 3-8% in nuclear compared 
to p-p collisions at the same c.m. energies, due to an overall net gluon antishadowing, although different regions of 
the top-quark differential distributions are depleted due to shadowing and EMC-effect corrections. The total cross 
sections for single-top, including the sum of t- and s-channels plus associated t W processes, are a factor of two 
(four) smaller than that for top-pair production at the LHC (FCC) and feature minimal nuclear modifications (±2% 
depending on the energy). 

After applying typical acceptance and efficiency cuts in the leptonic final-state, tt —> W~^b W~b —>■ b bii vv, 
one expects about 100 and 300 (anti)top-quarks per LHC-year and 5x10^ and 10^ per FCC-year at the nomi¬ 
nal luminosities in Pb-Pb and p-Pb collisions respectively. At the end of the LHC heavy-ion programme, with 
~ 10 nb“^, 1 pb“^ integrated Pb-Pb and p-Pb luminosities, about 2.5 thousand (fully-leptonic) f,f-quarks 
should have been measured individually by the CMS and ATLAS experiments. The number of visible single-top 
quarks produced in association with a W boson, in the similar tW W bW —?> bUvv final state, is lower by a 
factor of about 30 compared to tt production, due to the combination of lower cross sections, smaller reconstruction 
efficiencies, and only one top-quark per event. 

The proposed top-quark measurements at the LHC and FCC would not only constitute the first observation in 
nuclear collisions of the heaviest-known elementary particle, but would open up interesting novel physics opportu¬ 
nities such as constraints of nuclear parton densities in an unexplored kinematic range, studies of the dynamics of 
heavy-quark energy loss in the QGP, and colour-reconnection effects on the top-quark mass. We have, in particular, 
quantified the impact on the nuclear PDFs of the rapidity-differential distributions of the decay leptons from top- 
quark pairs, through the Hessian reweighting technique, finding that the data can be used to reduce the uncertainty 
on the Pb gluon density at high virtualities by up to 30% using the full LHC heavy-ion programme, and by about 
70% with just one FCC-year. 
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